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1 Introduction 

The theory of circle packings and, more generally, of circle patterns enjoys in recent years a fast 
development and a growing interest of specialists in complex analysis and discrete mathematics. 
This interest was initiated by Thurston's rediscovery of the Koebe-Andreev theorem about 
circle packing realizations of cell complexes of a prescribed combinatorics and by his idea about 



approximating the Riemann mapping by circle packings (see |T1, |RS[ ) . Since then many other 
remarkable facts about circle patterns were established, such as the discrete maximum principle 
and Schwarz's lemma and the discrete uniformization theorem |BS|. These and other results 



demonstrate surprisingly close analogy to the classical theory and allow one to talk about an 



emerging of the "discrete analytic function theory" [DS|, containing the classical theory of analytic 
functions as a small circles limit. 

Approximation problems naturally lead to infinite circle patterns for an analytic description of 
which it is advantageous to stick with fixed regular combinatorics. The most popular are hexagonal 
packings where each circle touches exactly six neighbors. The C°° convergence of these packings 



to the Riemann mapping was established in [HS|. Another interesting and elaborated class with 



similar approximation properties to be mentioned here are circle patterns with the combinatorics of 
the square grid introduced by Schramm . The square grid combinatorics of Schramm's patterns 
results in an analytic description which is closer to the Cauchy-Riemann equations of complex anal- 
ysis then the one of the packings with hexagonal combinatorics. Various other regular combinatorics 
also have similar properties Q]. 

Although computer experiments give convincing evidence for the existence of circle packing 
analogs of many standard holomorphic functions [ pS| ], the only circle packings that have been de- 
scribed explicitly are Doyle spirals | BDS[| (which are analogs of the exponential function) and con- 



formally symmetric packings ||BH|] (which are analogs of a quotient of Airy functions). Schramm's 
patterns are richer with explicit examples: discrete analogs of the functions exp(2;), erf (z), Airy |^ 
and z'^,log{z) [ |AB| ] are known. Moreover erf(z) is also an entire circle pattern]^ 

A natural question is: what property is responsible for this comparative richness of Schramm's 
patterns? Is it due to the packing - pattern or (hexagonal - square) combinatorics difference? Or 
maybe it is the integrability of Schramm's patterns which is crucial. Indeed, Schramm's square 



grid circle patterns in conformal setting are known to be described by an integrable system | BP2 
whereas for the packings it is still unknown^. 
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Doyle conjectured that the Doyle spirals are the only entire circle packings. This conjecture remains open. 
^It should be said that, generally, the subject of discrete integrable systems on l a ttices dif ferent from is 



underdeveloped at present. The list of relevant publications is almost exhausted by |Ad, KN, N 
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In the present paper we introduce and study hexagonal circle patterns with constant angles, 
which merge features of the two circle patterns discussed above. Our circle patterns have the 
combinatorics of the regular hexagonal lattice (i.e. of the packings) and intersection properties of 
Schramm's patterns. Moreover the latter are included as a special case into our class. An example 
of a circle patterns studied in this paper is shown in Figure |^. Each elementary hexagon of the 
honeycomb lattice corresponds to a circle, and each common vertex of two hexagons corresponds to 
an intersection point of the corresponding circles. In particular, each circle carries six intersection 
points with six neighboring circles and at each point there meet three circles. To each of the three 
types of edges of the regular hexagonal lattice (distinguished by their directions) we associate an 
angle < < ^r, n = 1, 2, 3 and set that the corresponding circles of the hexagonal circle pattern 
intersect at this angle. It is easy to see that a^'s are subject to the constraint ai + 02 + as = vr. 

We show that despite of the different combinatorics the properties and the description of the 
hexagonal circle patterns with constant angles are quite parallel to those of Schramm's circle pat- 
terns. In particular, the intersection points of the circles are described by a discrete equation of 
Toda type known to be integrable [Adj. In Section ^ we present a conformal (i.e. invariant with re- 



spect to Mobius transformations) description of the hexagonal circle patterns with constant angles 
and show that one can vary the angles a„, arbitrarily preserving the cross-ratios of the intersection 
points on circles, thus each circle pattern generates a two-parameter deformation family. Analytic 
reformulation of this fact provides us with a new integrable system possessing a Lax representation 
in 2 by 2 matrices on the regular hexagonal lattice. 

Conformally symmetric hexagonal circle patterns are introduced in Section ^. Those are defined 
as patterns with conformally symmetric flowers, i.e. each circle with its six neighbors is invariant 
under a Mobius involution (Mobius 180" rotation). A similar class of circle packings was inves- 



tigated in [BH]. Let us mention also that a different subclass of hexagonal circle patterns - with 
the multi-ratio property instead of the angle condition - was introduced and discussed in detail 
in [ BHS| . In particular it was shown that this class is also described by an integrable system. 



Conformally symmetric circle patterns comprise the intersection set of the two known integrable 
classes of hexagonal circle patterns: "with constant angles" and "with multi-ratio property". The 
corresponding equations are linearizable and can be easily solved. 

Further in Section |6| we establish a rather remarkable fact. It turns out that Doyle circle 
packings and analogue hexagonal circle patterns with constant angles are build out of the same 
circles. Moreover, given such a pattern one can arbitrarily vary the intersection angles an preserving 
the circle radii. 

Extending the intersection points of the circles by their centers one embeds hexagonal circle 
patterns with constant angles into an integrable system on the dual Kagome lattice (see Section ^) . 
Having included hexagonal circle patterns with constant angles into the framework of the theory of 
integrable systems, we get an opportunity of applying the immense machinery of the latter to study 
the properties of the former. This is illustrated in Section ^, where we introduce and study some 
isomonodromic solutions of our integrable systems on the dual Kagome lattice. The corresponding 
circle patterns are natural discrete versions of the analytic functions z*^ and logz. The results of 
Section |8| constitute an extension to the present, somewhat more intricate, situation of the similar 



construction for Schramm's circle patterns with the combinatorics of the square grid [BP2, AB| 
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2 Hexagonal circle patterns with constant angles 



The present paper deals with hexagonal circle patterns, i.e. circle patterns with the combinatorics 
of the regular hexagonal lattice (the honeycomb lattice). An example is shown in Figure |l[ Each 
elementary hexagon of the honeycomb lattice corresponds to a circle, and each common vertex 
of two hexagons corresponds to an intersection point of the corresponding circles. In particular, 
each circle carries six intersection points with six neighboring circles and at each intersection point 
exactly three circles meet. 

For analytic description of hexagonal circle patterns we introduce some convenient lattices and 




Figure 1: Doyle hexagonal isotropic circle pattern. 

variables. First of all we define the regular triangular lattice TC as the cell complex whose vertices 
are 

y(r£) = = + + ma;2 : A:,£,m G z|, where w = exp(27ri/3), (1) 

whose edges are all non-ordered pairs 

E{TC) = {[31,32] : 31,32 G V{TC), I31 -32I = l}, (2) 

and whose 2-cells are all regular triangles with the vertices in V{TC) and the edges in E{TC). We 
shall use triples {k,£,m) G as coordinates of the vertices. On the regular triangular lattice two 
such triples are equivalent and should be identified if they differ by the vector (n, n, n) with n G Z. 

The vertices of the regular triangular lattice correspond to centers and intersection points of 
hexagonal circle patterns. Associating one of the centers with the point k = £ = m = we obtain 
the regular hexagonal sublattice TCC with all 2-cells being the regular hexagons with the vertices in 

V{nC) = ^^ = k + iuj + mu? : k,i,meZ, k + i + m^O (mod 3)| , (3) 

and the edges in 

E{n£) = {[u,d2] ■■ 31,32 Gl^(W£), |3i-32| = l}. (4) 

The cells and the vertices correspond to circles and to the intersection points of the hexagonal circle 
patterns respectively. Natural labelling of the faces 

F{n£) = \} = k + euj + mu? : k,£,m£Z, k + £ + m = (mod 3)| 
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yields 

v{nc)\jF{nc) = v{TC). 



Definition 2.1 We say that a map w : V{7iC) ^ C defines a hexagonal circle pattern, if the 
following condition is satisfied: 

• Let 

Ik = l' + ^V{1-LC), /c = 1, 2, . . . , 6, where e = exp(7ri/3), 

he the vertices of any elementary hexagon inTCC with the center^' G F{7iC). Then the points 
w{ii),w{i2)-, ■ ■ ■ G C lie on a circle, and their circular order is just the listed one. We 

denote the circle through the points w{ii),w[i2)i ■ ■ ■ ■,w{iq) by C{^'), thus putting it into a 
correspondence with the center 3' of the elementary hexagon above. 

As a consequence of this condition, we see that if two elementary hexagons of TCC with the 
centers in 3', 3" G F{Ti.C) have a common edge [31,32] £ E{Ti.C), then the circles C{^') and C{^") 
intersect in the points ti^(3i), ^^(32)- Similarly, if three elementary hexagons of TCC with the centers 
in 3', 3", 3"' G Fine) meet in one point 30 G then the circles C{^'), C{f) and C{f) also 

have a common intersection point w{^o). 

Remark. We will consider also circle patterns defined not on the whole of TCC, but rather on 
some connected subgraph of the regular hexagonal lattice. 

To each pair of intersecting circles we associate the corresponding edge e G EiTCC) and denote 
by (l){z) the intersection angle of the circles, < (/> < 27r. The edges of EiTCC) can be decomposed 
into three classes 

E? = {t = [^',^"]GE{■HC)■.^'-f = ±l}, 

E^ = {t=[i',i"]eE{nC):i'-i" = ±u;], (5) 



E3 



^ = {z=b',f]eEinC):i'-f = ±u;'}. 



These three sets correspond to three possible directions of the edges of the regular hexagonal lattice. 

We shall study in this paper a subclass of hexagonal circle patterns intersecting at given angles 
defined globally on the whole lattice. 

Definition 2.2 We say that a map w : V(TCC) 1— > C defines a hexagonal circle pattern with 



constant angles if in addition to the condition of Definition \2.1\ the intersection angles of circles 
are constant within their class E^ , i.e. 

(j){e) = an Ve G n = 1,2,3. 

This angle condition implies in particular 

ai + 02 + "3 = TT. (6) 

We call the circle pattern isotropic if all its intersection angles are equal ai = 02 = «3 = ^r/S. 

The existence of hexagonal circle patterns with constant angles can be easily demonstrated via 
solving a suitable Cauchy problem. For example, one can start with the initial circles C{n{l — 
Lo)), C{n{uj - oj^)), C{n{uj^ - 1)), n G N. 
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Remark. In the special case ai = 02 = 7r/2,a3 = at each vertex one obtains two pairs 
of touching circles intersecting orthogonally. The hexagonal circle pattern becomes in this case 
a circle patterns of Schramm |Q with the combinatorics of the square grid. The generalization 
CKi + «2 = 71", «3 = was introduced in BP2 |. 

Figure |^ shows a nearly Schramm pattern. 



Figure 2: A nearly Schramm pattern. 



3 Point and radii descriptions 

In this paper three different analytic descriptions are used to investigate hexagonal circle patterns 
with constant angles. Obviously, these circle patterns can be characterized through the radii of the 
circles. On the other hand, they can be described through the coordinates of some natural points, 
such as the intersection points or centers of circles. Finally, note that the class of circle patterns 
with constant angles is invariant with respect to arbitrary fractional-linear transformations of the 
Riemann sphere C (Mobius transformations). Factorizing with respect to this group we naturally 
come to a conformal description of the hexagonal circle patterns with constant angles in Section ^ 
The basic unit of a hexagonal circle pattern is the flower, illustrated in Figure ^ and consisting 
of a center circle surrounded by six petals. The radius r of the central circle and the radii r„, n = 
1, . . . , 6 of the petals satisfy 

6 

arg JJ(r + e^""r„) (7) 

n=l 

where an is the angle between the circles with radii r and r„. Specifying this for the hexagonal 
circle patterns with constant angles and using (|6|) one obtains the following theorem. 

Theorem 3.1 The mapping r : F{7iC) is the radius function of a hexagonal circle pattern 

with constant angles 01,02, 03 iff it satisfies 

3 

argJ](l + e'""i?„)(e-'""+i?„+3) = 0, Rn = ^ (8) 

n=l 
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Figure 3: Circle flower. 



for every flower. 

Conjugating (|8|) and dividing it by the product Ri . . . Rq one observes that for every hexagonal 
circle pattern with constant angles there exists a dual one. 

Definition 3.2 Let r : F{7iC) be the radius function of a hexagonal circle pattern CP with 

constant angles. The hexagonal circle pattern CP* with the same constant angles and the radii 
function r* : F{TCC) — > M+ given by 

r* = I (9) 

is called dual to CP. 

For deriving the point description of the hexagonal circle patterns with constant angles it is 
convenient to extend the intersection points of the circles by their centers. Fix some point Pqo £ C. 
The reflections of this point in the circles of the pattern will be called conformal centers of the circles. 
In the particular case Poo = oo the conformal centers become the centers of the corresponding 
circles. We call an extension of a circle pattern by conformal centers a center extension. 

Mobius transformations play crucial role for the considerations in this paper. Recall that the 
cross-ratio of four points 

. N {Z2 - Zl){z4 - Zs) 

q[Zl,Z2, Z3, Z4) ■■= 7 TT r (10) 

(^3 - Z2){zi - Z4) 

is invariant with respect to these transformations. We start with a simple 

Lemma 3.3 Let Z2,Z4 be the intersection points and zi,Z3 the conformal centers of two circles 
intersecting with the angle a as in Figure The cross-ratio of these points is 

g(zi,Z2,23,-Z4) = (11) 

The claim is obvious for the Euclidean centers. These can be mapped to conformal centers by an 
appropriate Mobius transformation which preserves the cross-ratio. 



6 



Figure 4: Cross-ratio of an elementary quadrilateral. 



Thus a circle pattern provides a solution z : V{TC) ^ C to with the corresponding angles 
a. This solution is defined at the vertices of the lattice TC with quadrilateral sites comprised 
by the pairs of points as in Lemma |3.3| . For a hexagonal circle pattern with constant angles one 
can derive from (^) equations for the intersection points z : V{TCC) C and conformal centers 
z : V{TC \ nC) C. 

Theorem 3.4 Let z, zi,Z2, Z3, Z4, Z5, zq be conformal centers of a flower of a hexagonal circle pat- 
tern with constant angles 01,02,03, where an,n = 1,2,3 are the angles of pairs of circles corre- 
sponding to z,Zn and z,Zn+3- Define 81,82,63 through 

2an = 5n+2 - <5n+i, (mod 27r), n G {1, 2, 3} (mod 3). (12) 

Then z^ satisfy a discrete equation of Toda type on the hexagonal lattice 

j2Aj^— + —^]=Q, (13) 

where 

An = e*^"+2 - e^-^^+i n G {1, 2, 3} (mod 3). 

Let wi,W2,W3 be the neighboring intersection points to the point w of a hexagonal circle pattern 
and let Un be the angle between the circles intersecting at w,Wn- Then the following identity holds 




Figure 5: Intersection points and conformal centers of a Schramm circle pattern. 

In the special case of Schramm's patterns 03 = 0,ai + 02 = flowers contain only four petals and 
we arrive at the following 
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Theorem 3.5 The intersection points w^Wr-,Wu-,wuWd S C and the conformal centers z, Zr, Zu, zi, G 
C of the neighboring circles of a Schramm pattern (labelled as in Figure ^ satisfy the discrete equa- 
tion of Toda type on Tl? lattice 

1111 

+ = + , (15) 



W — Wr W — Wl W — Wu W — Wd 

^ + ^ = ^ + (16) 

Z — Zr Z — Zg Z — Zu Z — Zd 

The proofs of these two theorems are presented in Appendix ^ (in the case of complex a S C) . 

It should be noticed that both equations (|l6|) and ([l^) equations appeared in the theory of 
integrable equations in a totally different context |Su, |Ad|] . The geometric interpretation in the 
present paper is new. 

The sublattices of the centers and of the intersection points are dependent and one can be 
essentially uniquely reconstructed from the other. The corresponding formulas, which are natural 



to generalize for complex cross-ratios (IC), hold for both discrete equations of Toda type discussed 
above. We present these relations, which are of independent interest in the theory of discrete 
integrable systems, in Appendix ^ 

In Section |^ we will show that the hexagonal circle patterns with constant angles are described 
by an integrable system on a regular lattice closely related to the lattices introduced in Section |2|. 

4 Conformal description 

Let us now turn our attention to a conformal description of the hexagonal circle patterns. This 
description will be used in the construction of conformally symmetric circle patterns in the next 
section. We derive equations for the cross-ratios of the points of the hexagonal lattice that allow 
to reconstruct the lattice up to Mobius transformations. 

First we will investigate the relations of cross-ratios inside one hexagon of the hexagonal lattice 
as shown in Fig. |^. 




Figure 6: Cross-ratios in a hexagon. 



Lemma 4.1 Given a map w : V{7iC) — > C /ei 31, . . . ,3e be six points of a hexagon cyclically ordered 
(see Fig. To each edge of the hexagon, let us assign a cross-ratio Ti of successive points 

Wi = w{ii): 

Ti := q{wi,Wi+i,Wi+2,Wi-i), (i mod 6). (17) 

Then the equations 

Ti Tq T2 I-T2 ^ ' 

hold. 
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Proof. Let mj be the Mobius transformation that maps Wi-i,Wi, and lUj+i to 0, 1, and oo. Then 
Mi := m^^^rrii maps Wi-i,Wi, and Wj+i to Wi,Wi+i, and respectively and has the form 



1 -1 

1 -T,; 



Now {MaM^M^y^ = PM3M2M1 for some p gives the desired identities. □ 
Since every edge in E{7iC) belongs to two hexagons it carries two cross-ratios in general. 
We will now investigate the relation between them and show that in case of a hexagonal circle 
pattern with constant angles the two cross-ratios coincide. This way Lemma LI will furnish a map 
T : E{nC) C. 

Besides the cross-ratios of successive points in each hexagon (the Tj) we will need the cross-ratios 
of a point and its three neighbors: 

Definition 4.2 Let 31,32) o.nd 33 he the neighbors of ^ £ V{7iC) counterclockwise ordered and 
G ^i = 1,2,3. Any map w : V{7{C) C furnishes three cross-ratios in each point 
3 G V{nC): 

S^'^ := q{wi,Wi+i,w,Wi+2), i = l,2,3 ( mod 3). (19) 
They are linked by the modular transformation 

^r'^ = 1 - 4) ' 3) (20) 




Figure 7: Six points around one edge. 



Of course the two types of cross-ratios are not independent: When we look at an edge from 
E{Ti.C) with its four neighboring edges the six points form four cross-ratios (two T and two S). 
Three points fix the Mobius transformation and the other three can be calculated from only three 
cross-ratios. Therefore one expects one equation: 

Lemma 4.3 Given w : V(7iC) C and 31, . . . ,36 o-s shown in Fig. and the cross-ratios 

Ri := q{wi,W2,W4,W3), R2 := q{w5,W2,WQ,W3) (21) 



Sil! ■■= q{w3,W4,W2,WQ), S^'^ := q{w2,W5,W3,wi) (22) 



with [32,33] £ . Then the following identity holds: 



1-T,-R2- sS ■ ^ ^ 
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Proof. Use that Tj and i?j are linked by the relation 



Ri ■■= _ rj. = q{wi,Wi+i,Wi-i,Wi+2) (24) 

and insert the definition of the cross-ratio. □ 
Everything written so far holds for any hexagonal lattice, but in case of the hexagonal circle 
patterns with constant angles we can calculate all S in terms of the only: 

Sf^ := q{wi,w,+i,w,Wi+2) = (25) 

sm ai+2 

if [3«)3] ^ ^i- This can be verified easily by applying a Mobius transformation sending 3 to infinity. 
The three circles become straight lines forming a triangle with angles tt — . The are now the 
quotients of two of its edges. 

In particular we get S^]^ = S^^J along all edges [31,32] G E{TlC). Thus for our hexagonal circle 
patterns the above formula ( p5|) implies that there is only one T per edge since we get from equation 

dD 

Ti = T2. (26) 



So Lemma 4.1 defines in fact a map T : EiTiC) 



Note that we can get back the ai from say S^^^ in equation ( p5|) by the following formulas: 

-o(l) 1 , 0(1) 0(1) 1 , -a(l) 

Now we can formulate the main theorem that states that the equations which link the cross- 
ratios in a hexagon plus the constantness of S^^^ on the hexagonal lattice describe for real-negative 
T a hexagonal circle pattern with constant angles up to Mobius transformations: 

Theorem 4.4 

1. Given S" G C and a map w : V{nC) C for which Sf^ = S for all 3 € V{HC). Then the 
cross-ratios Ti of Lemma |^. j| define a map T : EiTLC) C and obey equations (p7 



2. Conversely given a solution T : E{H£) C of ^1^) then for each S G C there is up to 
Mobius transformations a unique map w : V{TCC) C having the map T as cross-ratios in 
the sense of Lemma ^TJ and having S^^^ = S for all 3 G ViTLC). 

If the solution T is real-negative the resulting w define a hexagonal circle pattern with constant 
angles. 

Proof. 1. is proven above. 

2. Starting with the points 11)2,0,01 ■u)i,i,0) and 11^1,0,1 (which fixes the Mobius transformation) we can 
calculate t«i,o,o using S^^qq. Now we have three points for each hexagon touching in 31,0,0- Using 
the T's we can determine all other points of them. Every new point has at least two determined 
neighbors so we can use the S^^^ to compute the third. Now we are able to compute all points of 
the new neighboring hexagons and so on. 

Since real-negative cross-ratios imply that the four points lie cyclically ordered on a circle, we 
have a circle pattern in the case of real-negative map T. But with (p5| ) the constantness of S^^^ 
implies that the intersection angles are constant too. □ 
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We have shown, that hexagonal circle patterns with constant angles come in one complex (or 
two real) parameter families: one can choose S £ C arbitrarily preserving T's. In Appendix ^ it 
is shown how this observation implies a Lax representation on the hexagonal lattice for the system 

dl)- 

Figure |2| shows a nearly Schramm pattern. One can obtain his description by taking combina- 
tions of Tj's and S^^^'s that stay finite in the limit Q3 — > 0. 

5 Conformally symmetric circle patterns 

The basic notion of conformal symmetry introduced in [ [BH[| for circle packings can be easily gen- 
eralized to circle patterns: Every elementary flower shall be invariant under the Mobius equivalent 
of a 180° rotation. 

Definition 5.1 1. An elementary flower of a hexagonal circle pattern with petals Ci is called 
conformally symmetric if there exists a Mobius involution sending Ci to Cj+3 (imod 3). 

2. A hexagonal circle pattern is called conformally symmetric if all of its elementary flowers are. 

For investigation of conformally symmetric patterns we will need the notion of the multi-ratio 
of six points: 

Zi- Z2Z3- Z4^Z5- Zq 

m{zi,Z2, Z3, Z4, Z5, zq) := . (28) 

Z2 - Z3Z4- Z5 ze - zi 



Hexagonal circle pattern with multi-ratio —1 are discussed in [ pHS (see p<!S| for further geometric 



interpretation of this quantity). It turns out that in the case of conformally symmetric hexagonal 
circle patterns the two known integrable classes "with constant angles" and "with multi-ratio —1" 
coincide. 

Proposition 5.2 1. An elementary flower of a hexagonal circle pattern is conformally symmet- 
ric if and only if the opposite intersection angles are equal and the six intersection points with 
the central circle have multi-ratio —1. 

2. A hexagonal circle pattern is conformally symmetric iff and only if it has constant intersection 
angles and for all circles the six intersection points have multi-ratio —1. 

Proof. First let us show that six points Zi have multi-ratio —1 if and only if there is a Mobuis 
involution sending Zi to Zi+3. 

If there is such Mobius transformation it is clear, that q{zi, Z2, z^, Z4) = q{z4, z^, zq, zi) and 

"1(^1, ^2, Z3, Z4,, Z5, Zq) = -q{zi, Z2, Z3, Z4,)/q{zi, Z5, Zq, Zi). (29) 

implies 

m{zi,Z2,Z3,Z4,Z5,ZQ) = -1. (30) 

Conversely let M be the Mobius transformation sending zi, Z2 and Z3 to Z4, z^ and zq. Then for := 
M{z4) equation (^) implies q{z4, Z5, zq, zi) = q{zi, Z2, Z3, Z4) = q{M{zi), M{z2), M{z3), M{z4)) = 
q{zi, Zq, z-t) and thus 3* = 31. The same computation yields 32 = (35) and 33 = M{iq). 

Now the first statement of the theorem is proven, since the intersection points and angles 
determine the petals completely. For the prof of the second statement the only thing left to show 
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is that all flowers being conformally symmetric implies that the three intersection angles per flower 
sum up to vr. So let us look at a flower around the circle C with petals Cj. Let the angle between 
C and Ci be (and we know that Oi = aj+s). Then the angle /3i = vr — ai — a2 is the angle 
between Ci and C2 and (S^ = n — — ai the one between C3 and C4. Since the flowers around C2 
and C3 are conformally symmetric too we now have two ways to compute the angle (32 between C2 
and C3. Namely 



TT 



^2 - "3 = P2 



TT 



(vr — as — ai) — (vr — ai — 02) 



which implies ai + 02 + «3 = tt- □ 
Using ( p9|) we see, that in the case of multi-ratio —1 the opposite cross-ratios T defined in 
section ^ must be equal: Tj = T^^i^^^ 3). Thus on EiTLC) the T's must be constant in the direction 
perpendicular to the edge they are associated with. For the three cross-ratios in a hexagon we get 
from (PD 



Ti + r2 + r3 - rir2T3 = 2. (31) 

Let us rewrite this equation by using the labelling shown in Fig. H. ^Ve denote by o^, hi and Cm 




Figure 8: Cross-ratios for conformally symmetric circle patterns 



the cross-ratios (^) associated with the edges of the families , E2 and E^ respectively. Note 
that for the labels in Fig. ^ 

k + (. + m = l 

holds. Written in terms of flfc, Cm equation ( |3l| ) is linear 



ak + bi + c„ 



(32) 



and can be solved explicitly: 



Lemma 5.3 The general solution to (3^) on E(7iC) is given by 

= bo + eA 
Cm = co + mA 



(33) 



with oo, bo, Co G C and A = 1 — qq — bo — cq. 
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Proof. Obviously ( pSf ) solves (^). On the other hand it is easy to show, that given the cross-ratios 
on three neighboring edges determine all other cross-ratios recursively. Therefore (33) is the only 
solution to (PI) □ 




Theorem 5.4 Conformally symetric circle hexagonal circle patterns are described as follows: Given 
an,bn,Cn by (35), choose S^^^ G C (or angles ai and then there is a conformally symmetric 



circle pattern with intersection angles ai,a2, and 03 given by formula 

The cross-ratio of four points can be viewed as a discretization of the Schwarzian derivative. 
In this sense conformally symmetric patterns correspond to maps with a linear Schwarzian. The 

shows that a symmetric solution of (5^) 



latter are quotients of two Airy functions [BH|. Figure 
is a good approximation of its smooth counterpart. Figure |2| also shows a conformally symmetric 
pattern with S^^^ = Wi and gq = bo = 1/3 — 0.29, and cq = 1 — 2aQ. 



6 Doyle circle patterns 

Doyle circle packings are described through their radii function. The elementary flower of a hexag- 
onal circle packing is a central circle with six touching petals (which in turn touch each other 
cyclically). Let the radius of the central circle be R and Ri, . . . ,Rq the radii of the petals. Doyle 
spirals are described through the constraint 

RkRk+3 = R^-, RkRk+2Rk+A = R^ ■ (34) 

There are two degrees of freedom for the whole packing — e. g. Ri/R and R2/R which are constant 
for all flowers. The next lemma claims that the same circles that form a Doyle packing build up a 
circle pattern with constant angles. 
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Lemma 6.1 The radii function of a Doyle packing (i. e. a solution to ^^)) solves the pattern radii 
equation Ij^) for any choice of the angles ai, 02, as = vr — ai — a2- 



Proof. Now insert identities (34) into (^). □ 

Definition 6.2 A hexagonal circle pattern with constant angles who's radii function obeys the 
constraint ^2) is called a Doyle pattern. 



Figure |T| shows a Doyle pattern. The fohowing lemma and theorem shows how the Doyle patterns 
fit into our conformal description: 

Lemma 6.3 Doyle patterns are conformally symmetric. 

Proof. We have to show that for Doyle pattern the multi-ratio of each hexagon is —1. One can 
assume that the circumfencing circle has radius 1 and center and that wi = 1. Then the other 
points are given by 

1 - Rji 



1 + i?,e*°0"'°'«3) 



Wj = Wi-l — 

where Rj are the radii of the petals. Inserting this into the definition of the multi-ratio implies the 
claim. □ 



Theorem 6.4 Doyle patterns and their Mobius transforms can be characterized in the following 
way: The corresponding solution to is constan^, i. e. + + cq = 1 and oq, 60) cq > 0. 

Proof. Any Doyle pattern gives rise to a constant solution of (^) since all elementary flowers of 
a Doyle pattern are similar. 

On the other hand one sees easily, that all < i?j < 1 can be realized. 

□ 



7 Lax representation and Dual patterns 

We start with a general construction of "integrable systems" on graphs which does not hang on 
the specific features of the lattice. This notion includes the following ingredients: 

• An oriented graph Q with the vertices V{Q) and the edges E[Q). 

• A loop group G[X\, whose elements are functions from C into some group G. The complex 
argument A of these functions is known in the theory of integrable systems as the "spectral 
parameter" . 

• A "wave function" ^ : V{Q) ^ G[A], defined on the vertices of Q. 

• A collection of "transition matrices" L : E{Q) ^ G[A] defined on the edges of Q. 

It is supposed that for any oriented edge e = {loutilin) G E{Q) the values of the wave functions in 
its ends are connected via 

^'(3,„,A) = L(e,A)4'(3o„t,A). (35) 

^Constant means here ak = ao,bk = 60, Cfc — co for all fc £ Z. 
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Therefore the fohowing zero curvature condition has to be satisfied. Consider any closed contour 
consisting of a finite number of edges of Q: 

fi = (31)32), e2 = (32,33), ep = (3p,3i)- 

Then 

L(ep,A)---L(e2,A)L(ei,A) = /. (36) 
In particular, for any edge e = (31,32) one has = (32,31) and 

L(e-\A) = (L(e,A))"\ (37) 

Actually, in applications the matrices L(e, A) depend also on a point of some set X (the "phase 
space" of an integrable system), so that some elements rc(e) G X are attached to the edges e of 
Q. In this case the discrete zero curvature condition (|3^) becomes equivalent to the collection of 
equations relating the fields x{zi), x{ep) attached to the edges of each closed contour. We 
say that this collection of equations admits a zero curvature representation. Such representation 
may be used to apply analytic methods for finding concrete solutions, transformations or conserved 
quantities. 




Figure 10: Quadrilateral lattice QC 



In this paper we will deal with zero curvature representations on the hexagonal lattice TiC 
and especially on a closely related to it a special quadrilateral lattice QC which is obtained from 
TiC by deleting from E{TC) the edges of the hexagonal lattice E{7iC) (i.e. corresponding to the 
intersection points). So defined lattice QC has quadrilateral cells, vertices V{QC) = V{TC) and 
the edges 

E{QC) = E{TC) \ E{nC) 
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as shown in Figure 10. The lattice dual to Q£ is known as Kagome lattice [Q. Similar to (|5|) there 
are tree types of edges in E[QC) distinguished by their directions: 

= {e = [3',3"]Gi^(Q/:): 3'-3" = ±^}, (38) 
E? = {t=[i',f]€E{QC):i'-f = ±co'}. 

There is a natural labelling {k,£,m) G of the vertices V{QC) which respects the lattice 
structure of QC Let us decompose V{Q/1) = Vq U Vi U V-i into three sublattices 

Vo = = k + luj + muj^ : k,i,m£Z, k + £ + m = 0^, 
Vi = ^^i = k + £uj + muj^ : k,£,meZ, k + £ + m = lY (39) 
V-i = = k + £uj + mu;^ : k,£,m£Z, k + £ + m = 

Note that this definition associates a unique triple {k, £, m) to each vertex. Neighboring vertices 
of QC, i.e. those connected by edges, are characterized by the property that their {k, £, m)-labels 
differ only in one component. Each vertex of Vq has six edges whereas the vertices of V±i have only 
three. 

The {k, £, m) labelling of the vertices of QC suggests to consider this lattice as the intersection of 
with the strip | /c + ^ + m |< 1. This description turns out to be useful especially for construction 

of discrete analogues of and logz in Section ^. We denote p = {k,£,m) G Tl? and three different 

types of edges of 1? by 

El = (p,p + (1,0,0)), £;2 = (p,p+ (0,1,0)), E3 = (p,P + (0,0,1)), pGZ^. 

The group G[A] we use in our construction is the twisted loop group over SL(2,C): 

G[A] = |l : SL(2,C) L(-A) = CT3L(A)ct3}, cig = diag(l, -1). (40) 

To each type En we associate a constant A„ G C. Let z^/^m^z'^im fields defined at vertices 
z, z* : ^ C and satisfying 

{Zin - Zout){z*n - Zl^t) = An (41) 

on ah the edges e = {pout,Pin) G En,n = 1,2,3. Here we denote Zout,in = z{pout,in)^ Kut,in = 
z*{pout,in)- We call the mapping z* dual to z. To each oriented edge e = {poutiPin) £ En we attach 
the following element of the group G[A]: 



lW(a) = (i-a2a„)-V2(^ 



f '^(^in Zqui) 

K4a 



(42) 

Note that this form as well as the condition (|4l| ) are independent of the orientation of the edge. 
Substituting (^ij) into ( ^2| ) one obtains L{X) in terms of the field z only. 

Dual fields z, z* can be characterized in their own terms. The condition that for any quadri- 
lateral of the dual lattice the oriented edges [z^^^, z^n) defined by ( |4l| ) sum up to zero imply the 
following statements. 
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Theorem 7.1 (i) Let z,z* : ^ C be dual fields, i.e. satisfying the duality condition (41)- Then 
three types of elementary quadrilaterals of 1? have the following cross-ratios: 



l{Zk/ 



Ai 
A3' 

Q[Zk/,m, Zk/,m+l, Zk/-l,m+l, Zk/-l,m) — -7—, 

( ^ - ^2 

^1 



(43) 



Same identities hold with z replaced by z* . 

(a) Given a solution z :W? ^ C to formulas determine (uniquely up to translation) 
a solution z* : ^ C of the same system (4-j)- 



It is obvious that the zero curvature condition (36) is fulfilled for every closed contour in if 
and only if it holds for all elementary quadrilaterals. It is easy to check that the transition matrices 
L*^") defined above satisfy the zero curvature condition. 

Theorem 7.2 Let z, z* : ^ C be a solution to and ei,e2,e3,e4 be consecutive positively 
oriented edges of an elementary quadrilateral of QC. Then the zero curvature condition 

L(e4,A)L(e3,A)L(e2,A)L(ei,A) = / 

holds with L(e, A) defined by 1^4^). Moreover let the elements 

LW(e,A) = (l-A2A„)-i/2(^ Y)' /5 = A„,n = l,2,3 

of G[X] be attached to oriented edges e = {pout,Pin) G En- Then the zero curvature condition on 1? 
is equivalent to existence of z, z* : 7^"^ ^ C such that the factorization 



/(e) = z{pin) - z{pout), 
holds. So defined z, z* satisfy (0, 



g(e) = z*{pin) - z*{pout) 



The zero curvature condition in Theorem is a generalization of the Lax pair found in [ N"C | 
for the discrete conformal mappings [BPl]. The zero curvature condition implies the existence of 
the wave function ^' : Z'^ G[X]. The last one can be used to restore the fields z and z*. There 
holds the following result having many analogs in the differential geometry described by integrable 
systems ("Sym formula", see, e.q., | BP2| ). 



Theorem 7.3 Let ^{p, A) be the solution of ( fgq j with the initial condition "^{p = 0, A) = /. Then 
the fields z, z* may be found as 



d^ 



k,£,m 



dX 



A=0 







^0,0,0 



Zk,e,m — 2:0,0,0 




(44) 



This simple observation turns out to be useful for analytic constructions of solutions in particular 
in Section ^. 

Interpreting the lattice QC as {{k, £, m) ^1? ■.\ /c + ^ + m|<l} one arrives at the following 



Corollary 7.4 Theorems \7.^ , \7.^ , |7. j| hold for the lattice QC, i.e. if in the statements one replaces 
1? by QC, p by 3, En by En and assumes k + i -\- m = in (4j). 
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Let us return to hexagonal circle patterns and explain their relation to the discrete integrable 
systems of this Section. To obtain the lattice QC we extend the intersection points of a hexagonal 
circle pattern by the conformal centers of the circles. The image of so defined mapping z : V(^QC) — > 
C consists of the intersection points z{y\ U V~\) and the conformal centers 2:(Vo). The edges 
connecting the points on the circles with their centers correspond to the edges of the quadrilateral 
lattice Q£. Whereas the angles a„ are associated to three types ^ n = 1, 2, 3 of the edges of the 
hexagonal lattice, constants (5„ defined in ( [T^ ) are associated to three types En,n = 1,2,3 of the 
edges of the quadrilateral lattice. Identifying them with 

A„ = e-'^- (45) 

in the matrices L^'^\X) above one obtains a zero curvature representation for hexagonal circle 
patterns with constant angles. 

Theorem 7.5 The intersection points and the conformal centers of the circles z : QC C of a 
hexagonal circle pattern with constant angles an,n = 1,2,3 satisfy the cross-ratio system on 
the lattice QC with An,n = 1,2,3 determined by (fl^, [^T 



The theorem follows from the identification of (^) with (^) using (|T^ ) and (|45|). 

The duality transformation introduced above for arbitrary mappings satisfying the cross-ratio 
equations preserves the class of such mappings coming from circle patterns with constant angles if 
one extends the intersection points by the Euclidean centers of the circles. The last one is nothing 
but the dual circle pattern of Definition 



Theorem 7.6 Let z : QC C be a hexagonal circle pattern CP with constant angles extended by 
the Euclidean centers of the circles. Then the dual circle pattern CP* together with the Euclidean 
centers of the circles is given by the dual mapping z* : QC C. 

Since A„, n = 1, 2, 3 are unitary (p5|), relation ^ follows directly from (]4l]). The patterns have the 



same intersection angles since the cross-ratio equation (43) for z and z* coincide. 



8 and log 2; patterns 

To construct hexagonal circle patterns analogs of holomorphic functions z'^ and log 2: we use the 
analytic point description presented in Section 0. Recall that we interpret the lattice QC of the 
intersection points and the centers of the circles as a subset \ k + £ + m \ < 1 oi and thus study 
the zero-curvature representation and ^'-function on Z^. 

A fundamental role in the presentation of this Section is played by a non- autonomous constraint 
for the solutions of the cross-ratio system (^) 

, 2 , ,7 r./7 s{zk+ll,m — Zkl,m){zkl,m — Zk-ll,m) , 

(>Zk,e,m + czk/,m + d= 2{k- ai) h 

Zk+l/,m — Zk-l/,m 

){Zk/ ,m ZjfP—lrn) /■An\ 

2{t - 02) 1- (46) 

Zk,e+l,m — Zk/-l,m 
„/ X {Zk,e,m+1 — Zk,e,m){Zk,e,m — Zk,i,m~l) 

2{m — 03) — — — — — . 

Zk/,m+l — Zk,e,m-1 

where 6, c, d, oi, 02, 03 G C are arbitrary. Note that the form of the constraint is invariant with 
respect to Mobius transformations. 
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Our presentation in this section consists of three parts. First, we explain the origin of the 
constraint ( ^6|) deriving it in the context of isomonodromic solutions of integrable systems. Then 
we show that it is compatible with the cross-ratio system (^3|), i.e. there exist non-trivial solutions 
of ( ^J4^ ). And finally we specify parameters of these solutions to obtain circle pattern analogs of 
holomorphic mappings z'^ and log z. 

We chose a different gauge of the transition matrices to simplify formulas. Let us orient the 
edges of 1? in the directions of increasing k ^ i ^ m. We conjugate L^^\X) of positively oriented 
edges with the matrix diag(l,A), and then multiply by (1 — A„A^)^/^ in order to get rid of the 
normalization of the determinant. Writing then for A^, we end up with the matrices 

(1 ^in ^out \ 

^ A„ ^ . (47) 

associated to the edge e = (pout^Pin) S En oriented in the direction of increasing of k -\- i + m. 
Each elementary quadrilateral of has two consecutive positively oriented pairs of edges Ci,t2 
and 63,64. The zero-curvature condition turns into 

Then the values of the wave function <I> in neighboring vertices are related by the formulas 

^k+i,e,m{fJ') = C^^\t,fi)^kAm{lJ'), e G El, 

< ^k,e+l,m{^^) = ^^'^\^,^^)^k,£,m{^^)^ eG^2, (48) 

{p)= C^'^\t,^i)^kAM, CG^3- 



We call a solution z : C of the equations (^3|) isomonodromic (cf. [^), if there exists a wave 

function $ : GL(2, €)[//] satisfying (H) and some linear differential equation in ji: 

J^^k/,m{fJ') = Ak,i,m{lJ')^k,i,mifJ'), (49) 

where Ak,e,mif^) are 2x2 matrices, meromorphic in /x, with the poles whose position and order do 
not depend on k, i, m. 

It turns out that the simplest non-trivial isomonodromic solutions satisfy the constraint (^). 
Indeed, since det C^^\^) vanishes at ^ = 1/A„ the logarithmic derivative of ^{fJ.) must be singular 
in these points. We assume that these singularities are as simple as possible, i.e. simple poles. 



Theorem 8.1 Let z : W? ^ C he an isomonodromic solution to with the matrix Ak/,m i'^ 

dT^ j of the form 

Ak,,,M = ^ + (50) 

^ n=l ^ - — 

with ix-independent matrices Ck/^mi ^ and normalized^ trace tr^o,o,o(/u) = 0. Then these 
®As explained in the proof of the theorem this normahzation can be achieved without loss of generality. 
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matrices have the following form: 

r - }i ( -^^k,e,m - c/2 bzl^^^^ + CZk/,m + d 

OZk/,m + c/2 



k,£,m 




b 


(1) 


k 


- a-i 


k,£,m 


Zk+l,£,m 




(2) 


e 


- 02 


k,£,m 


Zk,£+l,m 


^k,l—l,m 


(3) 


m 


- 03 


k,£,m 


Zk/,m+l 





1 Zk,£,m Zk~l,£,m ) 2 



oVJ; "o I Zk,£,m+\ ~ Zk,£,m {Zk,£,m+1 ~ Zk,£,m ,m Zk,t,m—l) \ I '^3 j 



and Zk/^m satisfies n4(\) - 

Conversely, any solution z : ^ C to the system ('4^,46) is isomonodromic with Ak,e,m{lJ-) 
given by the formulas above. 

The proofs of this theorem as well as of the next one are presented in Appendix |^. We prove 



Theorem by computing the compatibility conditions of (48) and ( |4g| ) and Theorem 8.2 by 
showing the solvability of a reasonably posed Cauchy problem. 

Theorem 8.2 For arbitrary b,c,d,ai,a2,a^ G C the constraint is compatible with the cross- 



ratio equations (4^); i-^- there are non-trivial solutions of ( [^J^ j. 
Further we deal with the special case of (^) where 

6 = ai = 02 = 03 = 0. 

If c 7^ one can always assume d = in ([46| ) by shifting z ^ z — —: 

„, {Zk+l,£ ,rn ^fc,£,m 

czk,e,m = 2 A; h 

Zk+l,e,m ~ Zk-l,£,rri 
r,/i\Zk,£+l,m ~ Zk,£,m ,m Zk/-l,m) _|_ ^^^^ 

Zk,£+l,m — Zk,£-l,m 
r, {Zk,£,m+1 — Zk,e,m){zk,e,ni — Zk,£,m-l) 

2m — . 

Zk,£,m+l Zk^£^m—1 

To define a circle pattern analogue of z'^ it is natural to restrict the mapping to the following two 
subsets of Z^: 

Q = {{k,£,m) eZ^ \ k>0,£>0,m <0}, H = {{k,e,m) eZ^ \ m <0} CZ^. 



Corollary 8.3 The solution z : Q ^ C of the system satisfying the constraint ( pl\ ) is uniquely 
determined by its values 

zi,o,o, ZQ^ifi, 2;o,o,-i- (52) 

Proof. Using the constraint one determines the values along the coordinate lines Znflfi, ZQ^nfl, -^co.-n 
Vn E N. Then all other Zk/^rm k,£, —nT- S N in consecutive order are determined through the cross- 
ratios (|^). Computations using different cross-ratios give the same result due to the following 
lemma about the eighth point. 
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Lemma 8.4 Let H{pi + l/2,p2 + \/2,p^ — 1/2), p = {pi,P2,P3) £ be the elementary hexahedron 
(lying in C) with the vertices ■Zp+(fc,£,-m)) k,i,m £ {0,1} and let the cross-ratios of the opposite 
cites of H{pi + 1/2, p2 + 1/2, — 1/2) be equal 



<l{Zp, Zp+(o,o,-l)> 2:p+(o,i,-i), 2p+(o,l,0)) = 9(^p+(l,0,0)i ^p+{l,0,-l)i ^p+{l,l,0)) 
liZp, 2;p+(l,o,0)> ^p+(l,0,-l)) ^p+(0,0,-l)) = 9(^p+(0,l,0); ^p+(l, 1,0)) ^p+(l, 1,-1)' ^P+(0,1-1)) 
<l{Zp, ^p+(0,l,0)> ^p+(l,l,0)i ^p+{l,0,0) ) = 9(^p+(0,0,-l)> ^p+(0,l,-l)> ^p+(l, 1,-1) > ^p+(l, 0,-1)) 



Q2 
93 



with 91^293 = 1- Then all the vertices of the hexahedron are uniquely determined through four given 

Zp, Zpj^{lfifl), ^p+(0,l,0); 2p+(0,0,-l)- 



The relation of solutions of (43,51) to circle patterns is established in the following 
Theorem 8.5 The solution z : Q ^ C of the system (J^,51) with the initial data 

2:1,0,0 = 1, ^0,1,0 = e*^, 2^0,0,-1 = e*'^ (53) 

and unitary cross-ratios qn = e~'^*"" determines a circle pattern. For all {k,£,m) G Q with even 
k + I + m the points Zk±i/^m, Zk/±i,m, Zk,i,m±i He on a circle with the center Zk,i,m- 

Proof. We say that a quadrilateral is of the kite form if it has two pairs of equal adjacent edges, 
(and that a hexahedron is of the kite form if all its cites are of the kite form) . These quadrilaterals 
have unitary cross-ratios with the argument equal to the angle between the edges (see Lemma p.3| ). 
Our proof of the theorem is based on two simple observations: 

(i) If a quadrilateral has a pair of adjacent edges of equal length and unitary cross-ratio then it 
is of the kite form, 

(ii) If the cross-ratio of a quadrilateral is equal to q{zi, Z2, Z3, z^) = e~^*" where a is the angle 
between the edges {zi,Z2) and (22,-^3) (as in Figure Q) then the quadrilateral is of the kite 
form. 



Applying these observation to the elementary hexahedron in Lemma 8.4 one obtains 



Lemma 8.6 Let three adjacent edges of the elementary hexahedron in Lemma \8.4 have equal length 
and the cross-ratios of the cites be unitary qn = e^^*"". Then the hexahedron is of the kite form. 



The constraint (51) with i = m = implies by induction 

I ^2n+l,0,0 — ^2n,0,0 | = | 2;2n,0,0 — ^2n,-l,0,0 I) 

and all the points 2;n,o,o lie on the real axis. The same equidistance and straight line properties hold 
true for i and m axes. Also by induction one shows that all elementary hexahedra of the mapping 
z : Q ^ C are of the kite form. Indeed, due to Lemma the initial conditions (|5^ ) imply that 
the hexahedron ff(l/2, 1/2, — 1/2) is of the kite form. Since the axis vertices lie on the straight 
lines, we get, for example, that the angle between the edges (zi,i,0) -2i,o,o) and (zi,o,0; -^2,0,0) is 03- 
Now applying the observation (ii) to the cites of H{3/2, 1/2, —1/2) we see that 

I ^2,0,0 — ^1,0,0 1 = 1 ^2,0,0 — ^2,1,0 1 = 1 ^2,0,0 — ^2,0,-1 | ■ 

Lemma p.6| implies H{3/2, 1/2, —1/2) is of the kite form. Proceeding further this way and control- 
ling the lengths of the edges meeting at z^^i^^ with even k -\- i -\- m and the angles at z/^^i^^n with 
odd k -\- i -\- m one proves the statement for all hexahedra. 
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Denote the vertices of the hexagonal grid by (see Figure 10) 

Qh = {{k,e,m)£Q:\k + i + m\<l}, 
Hh = {{k,i,ra) <^ H ■.\k + i + m\<l]. 

The half-plane Hh consists of the sector Qh and its two images under the rotations with the angles 
±7r/3. 

Corollary 8.7 The mapping z : Hh C given by (J^5l) with unitary cross-ratios qn = e~^*"" 
and unitary initial data 

I z±iflfi 1=1 2;o,±i,o 1=1 ^0,0,0 1= 1 

is a hexagonal circle pattern with constant angles an extended by the Euclidean centers of the circles. 
The centers of the corresponding circles are the images of the points with k + i + m = 0. 

The circle patterns determined by most of the initial data /3, 7 S M are quite irregular. But for a 
special choice of these parameters one obtains a regular circle pattern which we call the hexagonal 
circle pattern z'^ motivated by the asymptotic of the constraint (^) as k,£,m — > 00. 

Definition 8.8 The hexagonal circle pattern z^, < c < 2 (extended by the Euclidean centers of 
the circles) is the solution z : Hh C of ^1,51) with 

e^'"", n (mod 3), 



^'1+1 2m„ 



and the initial conditions 

21,0,0 = 1, ^0,1,0 = e ^ 2:0,0,-1 = e % 2_i,o,o = e , 20,-1,0 = e ^ (54) 



Motivated by our computer experiments (see in particular Figure 11) and the corresponding re- 
sult for Schramm's circle patterns with the combinatorics of the square grid | [AB| , Q we con- 
jecture that the hexagonal z'^ is embedded, i.e. the interiors of all elementary quadrilaterals 
(2^(31), 2^(32), ^(33), ^(34)), I 3i+i -3i 1= 1;* (mod 4) are disjoint. 

In the isotropic case when all the intersection angles are the same 

TT 

"1 = "2 = "3 = 2 • 

one has 

Ai = wAa = w^Ag, (55) 

and due to the symmetry it is enough to restrict the mapping to Qh- The initial conditions (U) 
become 

^1,0,0 = 1, ^0,1,0 = e2--/3, 20,0,-1 = (56) 

TT TT TT 

Like for the smooth z^ the images of the coordinate axes arg % = — , — , — are the axes with the 

vrc vrc vrc 
arguments ^^'^ ~^ respectively: 

vrc TTC TTC 

arg2„o,-n = arg2o,„,-n = -^7, arg2„,„,m = — , m = -2n; -2n± 1, Vn € N. 
2 6 
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Figure 11: Non-isotropic and isotropic circle patterns z^l"^ . 




Figure 12: Circle pattern z^l^ . 
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6 

For c = -, q £ N the circle pattern z'^{Qh) together with its additional q — I copies rotated 

by the angles 2TTn/q, n = 1, . . . ,q — 1 comprise a circle pattern covering the whole complex plane. 
This pattern is hexagonal at all circles beside the central one which has q neighboring circles. The 
circle pattern 

^6/5 jg 

shown in Figure O. 



Definition 8.8 is given for < c < 2. For c = 2 the radii of the circles intersecting the central 
circle of the pattern become infinite. To get finite radii the central circle should degenerate to a 
point 

2^0,0,0 = ^1,0,0 = ^0,0,-1 = ^0,1,0 = 0. (57) 

The mapping z : Q — > C is then uniquely determined by Z2,o,0) zo^2,o, 2^0,0,-2; -2i,i,0) ^1,0,-1, 2:0,1,-1, the 
values of which can be derived taking the limit c — > 2 — 0. Indeed, consider the quader m = with 
the cross-ratios of all elementary quadrilaterals equal to q = e"^*"^. Choosing 2^1,0,0 = 20,1,0 = ee*'^'^ 

as above we get zono = • Normalization 

c = 2 - 2e 

yields 

22,0,0 = 1, 20,2,0 = e'"''. (58) 

CO; 

Since the angles of the triangle with the vertices 20,0,0,21,0,0,21,1,0 are — respectively, 

one obtains 

p - sin(ca/2) 

21,1,0 — e + /leC , -Ke — e — r— . 

sin(aej 

In the limit e ^ we have 

Si'^" ia /CAN 

21,1,0 = e'". (59) 

a 

Observe that in our previous notations a = ai + 02 = tt — 03. Finally, the same computations for 
the sectors {£ = 0,k > 0,m < 0} and {k = 0, 1 > 0,m < 0} provide us with the following data: 

22,0,0 = 1, 20,0,-2 = e2*"2^ zo,2,o = e2^("i+°2)^ 

21,0,-1 = ^e-, .0,1,-1 = ^e^(— .1,1,0 = !i^^(^^^e^(".+-). (60) 
a2 Oil ai + a2 

In the same way the initial data for other two sectors of Hh are specified. The hexagonal circle 
pattern with these initial data and c = 2 is an analog of the holomorphic mapping . 
The duality transformations preserves the class of circle patterns we defined. 

Theorem 8.9 

{z^)* = z^\ c* = 2-c, 

where (z^)* is the hexagonal circle pattern dual to the circle pattern z*^ and normalized to vanish at 
the origin (z'^)*(3 = 0) = 0. 

Proof. Let us consider the mapping on the whole Q. It is easy to see that on the axes the duality 
transformation ( ]4l| ) preserves the form of the constraint with c being replaced by c* = 2 — c. Then 
the constraint with c* holds for all points of Q due to the compatibility in Theorem ^.2| . Restriction 
io\k + i + m\ <l implies the claim. 
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The smooth Umit of the duahty transformation of the hexagonal patterns is the fohowing 
transformation f ^ f* of holomorphic functions: 

{riz))' ^ 



f'iz) 



The dual of f{z) = is, up to constant, P{z) = logz. Motivated by this observation, we define 
the hexagonal circle pattern log z as the dual to the circle pattern z'^ : 



logz:= {z^y. 



2k h 2t h 

Zk+l,e,ni ~ Zk-l,e,m Zk/+l^m ~ Zk/-l,m 



The corresponding constraint (I 



„ {Zkl,m+1 — Zkl,m){zk,e,m — Zkl,m-l) ^ f„.^ 

Zm = 1. (blj 

Zk/,m+l — Zk,e,m-1 



can be derived as a limit c — > +0 (see [AB| for this limit in the square grid case). The initial data 
for logz are dual to the ones of z^. In our model case of the quader m = factorizing q = A1/A2 



with Ai = 1/2, A2 = e2*"/2 one arrives at the following data dual to (|^J5|,|^): 

1 

^0,0,0 — 00, ^^1,0,0 — 0) ^2,0,0 — 2' 

1 a 

zo,i,o = ia, zo,2,o = 7: + ^1,1,0 = TT-- (62) 

2 2 sm a 

In the isotropic case the circle patterns are more symmetric and can be described as mappings of 
Qh. 

Definition 8.10 The isotropic hexagonal circle patterns z^ and logz are the mappings Qh C 
with the cross-ratios of all elementary quadrilaterals equa^ to e~^'^*/^. The mapping z^ is deter- 
mined by the constraint ( fji^ with c = 2 and the initial data 

zq,o,o = 00, zififi = 0, zofi-i = 0, 20,1,0 = 0, 
^2,0,0 = 1, ^0,0,-2 = e^-/^ zo,2,o = 

zi,o -1 = — — e'^*/^ 20,1,-1 = — 2^1,1,0 = -— e^'^*/^ 
zir ZTT 47r 



For log z the corresponding constraint is (6J_) and the initial data are 



vr 27r 

^0,0,0 = 00, 21,0,0 = 0, 20,0,-1 = gi, 20,1,0 = —^, 

1 1 vr. 1 27r. 

22,0,0 — 2' ^0,0,-2 ~ 2 "3*' ~ 2^ "3"^' 

21,0,-1 = + ^), 20,1,-1 = ^(-^ + 3^), 21,1,0 = ^(-^ + ^)■ 

The isotropic hexagonal circle patterns z"^ and log 2 are shown in Figure |l^. 
Starting with 2^, c G (0,2] one can easily define z^ for arbitrary c by applying some simple 
transformations of hexagonal circle patterns. The construction here is the same as for Schramm's 



^The first point in the cross-ratio is a circle center and the quadrilaterals are positively oriented. 
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Figure 13: Isotropic circle patterns and logz. 



patterns (see Section 6 of [ |AE| ] for details) . Applying the inversion of the complex plane z i— > l/z 
to the circle pattern z^, c G (0, 2] one obtains a circle pattern satisfying the constraint with — c. It 
is natural to call it the hexagonal circle pattern z~'^, c G (0, 2]. Constructing the dual circle pattern 
we arrive at a natural definition of z'^'^'^. Intertwining the inversion and the dualization described 
above, one constructs circle patterns z^ for any c. 

In particular, inverting and then dualizing z = k + iuj + muj"^ with Ai = —3, A2 = — 3w^, A3 = 
—?>ui we obtain the circle pattern corresponding to z^: 

Zk/,m = {k + luj + muj^)^ — (A; + £ + m). 

Note that this is the central extension corresponding to Poo = 0. The points with even k + 1 + m 



can be replaced by the Euclidean centers of the circles. As it is shown in Section 6 of |AB] the 
replacement Poo = by Poo = 00 preserves the constraint (^l]) . 

9 Concluding remarks 

We restricted ourself to the analysis of geometric and algebraic properties of hexagonal circle 
patterns leaving the approximation problem beyond the scope of this paper. The convergence of 
circle patterns with the combinatorics of the square grid to the Riemann mapping is proven by 
Schramm in |^. We expect that his result can be extended to the hexagonal circle patterns defined 
in this paper. 

The entire circle pattern erfz also found in remains to be rather mysterious. We were 
unable to find its analogue in the hexagonal case, and thus, have no counter-examples to the Doyle 
conjecture for hexagonal circle patterns with constant angles. It seems that Schramm's erfz is a 
feature of the square grid combinatorics. 

The construction of the hexagonal circle pattern analogs of z^ and log z in Section 8| was based 
on the extension of the corresponding integrable system to the lattice 1? . Theorem 8^ claims that 



in this way one obtains a circle pattern labelled by three independent indices /c, ^, m. Fixing one 
of the indices, say m = niQ, one obtains a Schramm's circle pattern with the combinatorics of the 
square grid. In the same way the restriction of this "three-dimensional" pattern to the sublattice 
I k + £ + m — uq |< 1 with some fixed no G Z yields a hexagonal circle pattern. In Section |8| 
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we have defined circle patterns z'^ on the sublattices ttt-q = and ng = 0. In the same way the 
sublattices mo 7^ and no 7^ can be interpreted as the circle pattern analogs of the analytic 
function + a)*^, a 7^ with the square and hexagonal grids combinatorics respectively. 

It is unknown, whether the theory of integrable systems can be applied to hexagonal circle 
packings. As we mentioned already in the Introduction, the underdevelopment of the theory of 
integrable systems on the lattices different from may be a reason for this. We hope that 
integrable systems on the hexagonal and related to it lattices appeared in this paper and in [BHS 
as well as the the experience of integration of hexagonal circle patterns in general will be helpful 
for the progress in investigation of hexagonal circle packings. 



A Appendix. A Lax representation for the conformal description. 

We will now give a Lax representation for equations (|^). For 3 G Vi^C) let [31,3] G and be 
the Mobius transformation, that sends (31,32,33) to (0, 1, 00). For e = [3,3] G Ef set Li = m^om"^. 
The Li depend on 5^^^ and only. They are given in equation (|63|). Note that we do not need 
to orient the edges since LiL~^ is the identity if we normalize det-Lj = 1. The claim is that the 
closing condition when multiplying the Lj around one hexagon is equivalent to equation (|l8|): 



Theorem A.l Attach to each edge of Ei a matrix Li{T,S) of the form: 

Li{T,S) = 

L2{T,S) = { ^ r^a 1 (63) 

then the zero curvature condition for each hexagon in F{TCC) 

Li(T4, S)L2{n, S)L3{Te, S) = pL^in, 5)L2(T2, S)Li{Ti,S) (64) 
for all S is equivalent to ^T^). 

Proof. Straight forward calculations. □ 




Equations (63) and (|64|) are a Lax representation for equation (LS) with the spectral parameter S. 



B Appendix. Discrete equations of Toda type and cross-ratios 

In the following we will discuss briefly the connection between discrete equations of Toda type 
and cross-ratio equations for the square grid and the dual Kagome lattice. This interrelationship 
holds in a more general setting, namely for discrete Toda systems on graphs. This situation will 
be considered in a subsequent publication. 

We will start with the square grid. Let us decompose the lattice into two sublattices 

Zl = {{m, n) G (m + n mod 2) = k}, k = 0, 1. 



Theorem B.l Let z :1? ^ ^ be a solution to the cross-ratio equation 
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Then restricted to the sublattices Zq or 'L\ it satisfies the discrete equation of Toda type (|i 
1111 



+ 



^m+l,n+l 



^m— l,ra— 1 



+ 



^m+l,n—l ^m,n ^m— l,n+l 



(66) 



Conversely, given a solution z : Zq — > C o/ equation (66) and arbitrary q,w £ C there exists unique 
extension of z to the whole lattice z : I? ^ C with zi^ = w satisfying the cross-ratio condition 
(61). Moreover, restricted to the other sublattice the so defined mapping z : Zif ^ C satisfies 



the same discrete equation of Toda type (66). 



Proof. For any four complex numbers the cross-ratio equation can be written in a simple fraction 
form: 



q{ui,U2,U3,U4) =q<^ 



1 



1 



+ - 

U2 — Ui U2 — Us U2 — U4 



0. 



(67) 



To distinguish between the two sublattices, we will denote the points belonging to Zq with Zi while 



Z4 



W3 



Z3 



W4 









z 



V02 



Wi 



Z2 




Figure 14: Nine points of and 13 points of ViTC). 



Wi Are the points associated with Z^. Now (67) gives four equations for a point z and its eight 
neighbors as shown in Fig |l^: 

(9-1) 
- 1) 



(9-1) 



(9 



-1 



1 


z 






1 


z 


- Z2 




1 


z 


- Z3 




1 


z 


- Z4^ 



' Z — W4 Z — Wl 



-1 



1 



1 



Z — Wl Z — W2 
1 1 



(68) 



' Z — W2 Z — W3 



-1 



1 



1 



Z — W3 Z — W4 



Multiplying the second and fourth equation with q and taking the sum over all four equations (|6^) 



yields equation (6€). Conversely given wi, the first three equations y^) determine W2,W3, and W4 
uniquely. The so determined W3, and W4 satify the fourth equation ( p8| ) (for any choice of wi) if 
and only if ( |66|) holds. This proves the second part of the theorem. Obviously interchanging w and 
z implies the claim for the other sublattice. □ 



Now we pass to the dual Kagome lattice shown in Fig. 10, where a similar relation holds for the 
discrete equation of Toda type on the hexagonal lattice. However, the symmetry between the 
sublattices is lost in this case. 

Here we decompose V{TC) into V{nC) and V{T£) \ V{nC) ^ F{nC). 
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Theorem B.2 Given qi,q2,Q3 G C with qiq2Q3 = 1 md a solution z : V{TC) ^ C of the cross- 
ratio equations 

q{zk,e,m-i, Zk/,m, Zk,e+i,m, Zk,e+i,m-i) = qi 

,mi Zk/,mi Zk,e,ni-1: -2fe+l/,m-l ) = q2 (69) 
q{Zk,e-l,m, Zk,e,m, Zk+l,e,m, Zk+l,e-l,m) = ^3 

then 

1. restricted to F(TIC) the solution z satisfies the discrete equation of Toda type ( jl^ on the 
hexagonal lattice 

Y^Ak(^— + —^)=0 (70) 

^ \Z-Zk Z-Zk+3j 

2. restricted to V{TCC) the solution z satisfies 

k=l 

where Ai = Ajj_2 — Aj+i with Aj defined through qi = 

Conversely given (?i,(?2)93 £ C to^/i (?i(/2</3 = 1, a solution z to (%) on FiTiC) and 2^0,0,0; ^/^ere is 
a unique extension z : V{TC) C satisfying (6§). 



Given qi,q2,q3 G C with qiq2q3 = 1, a solution z to (71) on V(7iC) and -Zi,o,0; there is a unique 
extension z : V{TC) C satisfying Ml 



Proof. First ( pT] ) shows immediately that (^Tj) is equivalent to having constant 5 as described in 
Section |5[ 

Again we will distinguish the sublattices notationally by denoting points associated with ele- 
ments of V{Ti.C) with Wi and points from FiTiC) with Zi. If the cross-ratios and neighboring points 
of a z G F(TIC) are labelled as shown in Fig. |l^, equation ( |67| ) gives 6 equations: 

1 1 1 

{qi - 1) = qi- 



z — z\ z — wq z — Wi 
(91 - 1) = 91 

2 — ^4 2^ — ^)3 Z — 
(92 - 1) = 92 

2-^2 - z - W2 (72) 

(92 - 1) = 92 

z — z^ z — W4, z — 

(93 - 1) = 93 

Z - 2:3 Z -^W2 Z -^W3 

(93 - 1) = 93 

Z — Zq z — W5 z — Wq 

To proof the first statement we take a linear combination of the equations (|7^) . Namely a the first 
two plus b times the second two plus c times the third two. It is easy to see that there is a choice 
for a, b, and c that makes the right hand side vanish if and only if 51^293 = 1- If 9i = ; 92 = 
and 93 = ^ choose a = A2, 6 = A3, and c = Ai. The remaining equation is (70). 

To proof the third statement we note that given wi, we can compute W2 through wq from 
equations (72) 2. to 6. but the "closing condition" equation ( |72[ ) 1. is then equivalent to ([70|). 

The proofs of the second and fourth statements are literally the same if we choose Zi = ^^+3 
and Wi = Wi^3. □ 
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C Appendix. Proofs of Theorems of Section |^ 
Proof of Theorem |8.1| 

Let ^k/,ni{fJ') be a solution to (^,46) with some ^u-independent matrices Ck/,m, b'^^ ^. The deter- 
minant identity 

det <^>kAM = (1 - /«Ai)^(l - fiA^Yil - fiA^r det $o,o,o(/u) 

impUes 

... k £ m 

/i-l/Ai /i-l/A2 /i-l/Ag 

with a(//) independent of fc, ^, m. Without loss of generality one can assume a{n) = 0, i.e. 

tr b'i^I =k,tr B^i^] = I, tr B^'] 

This can be achieved by the change $ i— > exp(— 1/2 J a(fj,)dfi)^. 
The compatibility conditions of (Esl) and (p9|) read 



(73) 



dfi 


— Ak+i,e,m^^^^ 


- J^^'^Ak/, 


d/j, 




- C^^^Ak,i, 


dfi 


= Ak/,rn+l^^'^^ 





The principal parts of these equations in ^ = 1/Ai imply 

^kll,l,m [^}_ 1 j = ( J_ I \ ^klm^ h = Zk+l,e,m " Zk,t 



/2A 




mi 



The solution with tr B^^j ^ = k is 

p(l) _ k — ai / Zk+l^e^m — Zk/,m {Zk+l,£,m — Zk,e,m){Zk,e,m — Zk-l,e,m) \ I Qi r 

k Im — \i /o 

Zk+l/,m — Zk-l/,m \ ^ Zk/^m — ^fc-l/,m. / ^ 

l_ 1 (2) (3) 

The same computation yields the formulas of Theorem 3.1 for Bj^'^^ and -B^/^- 

To derive a formula for the coefficient Ck,t^m let us compare ^k,t,m{lA with the solution ^k,t,m{^) 
in Theorem 7.3, more exactly with its extension to the lattice Z'^: 



*fc+l,£,m(A) = L^^\\)^k/,^ 
^k/,m+l{\) = L^^\\)^k/,, 
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normalized by ^&o,o,o('^) = I- Here the matrices L*^") are given by (42) 
Consider 
with 

h{\) = (1 - - a2A2)^/2(1 - A^Aa)"/^. 

So defined is a function of /i. Since it satisfies the same difference equations ( p8[ ) as $(/i) and is 
normahzed by 

^k/At^ = 0) = /, 

we have 

^fcAm(Ai) = ^fcAm(M)$o,o,o(M)- (74) 
Moreover ^(/i) is holomorphic in ^ = and due to Theorem equal in this point to 

^k,e,m{fJ' ~ ~ Q l)' ^ ~ ^'^'^'"^ ~ ^0:0,0- 

Taking the logarithmic derivative of ( [7^ ) with respect to n 



and computing its singularity at /i = we get 



_. 1 Z\ ( 1 -Z 

'^k,t,m - 1 1 / I 1 



The formula for Ck^i^m in Theorem 8.1 is the general solution to this equation. 

Conversely, by direct computation one can check that the compatibility conditions ( |73|) with 
■A.k,i,m computed above are equivalent to (^6|). 



Proof of Theorem 8.2 



Proof. First note that one can assume 6 = by applying a suitable Mobius transformation. 
Next by translating z one can make -Zfc,i,m = for arbitrary fixed k,l,m £ Z (this will change 
d however). Finally we can assume d = or d = 1 since we can scale z. Now one can show 
the compatibility by using a computer algebra system like MATHEMATICA as follows: given the 
points Zk,i,m, Zk+i,i,m, Zk,i±i,m, Zk,i,m±i One can compute Zk-i,i,m using the constraint (||). With 
the cross-ratio equations one can now calculate all points necessary to apply the constraint for 
calculating Zk+2,i,m and Zk,i,m+2 (besides Zk-i,i,m-i one will need all points Zk±i,i±i,m±i) ■ Once 
again using the cross-ratios one calculates all points necessary to calculate -Zfc+i,/,m+2 with the 
constraint. Finally one can check, that the cross-ratio qizk,i,m+i, Zk+i,i,m+i, Zk+i,i,m+2, Zk,i,m+2) 
is correct. Since the three directions are equivalent and since all initial data was arbitrary (i.e. 
symbolic) this suffices to show the compatibility. □ 
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